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This  report  presents  an  ar.al7ti.cal  method  for  determining  the  flow 
characteristics  of  a  two-phase  flow  of  liquid  drops  in  a  gas  stream  as 
the  mixture  expands  through  an  annular  converging-diverging  nozzle. 

The  subject  analysis  car.  bt  utilized  to  predict  the  liquid  velocity,  gas 
velocity,  static  pressure,  .ini  droplet  diameter  as  a  function  of  axial 
distance  along  the  nozzle  for  a  two-phase  flow  that  contains  approxi¬ 
mately  ten  times  as  much  liauid  as  gas  by  weight. 

The  analysis  was  programmed  for  solution  on  an  RPC  i*000  Digital 
Computer.  Two  nozzle  configurations  were  investigated.  Roth  nozzles 
hed  the  same  converging  angle  of  20°,  throat  radius  of  one  inch,  inlet 
area  of  3-237  ir..B,  throat  area  of  0.267  in.^  and  exit  ar^a  of  2.150  in.* 
Oie  nozzle  had  a  total  diverging  angle  of  7°  and  the  other  had  s  total 
diverging  angle  of  21°, 

The  liquid  and  gas  utilized  in  this  investigation  were  water  and 
air  respective’-'.  Flow  rates  of  between  8  and  11  lb/sec  of  water  and 
between  1.0  and  1.3  lb/sec  of  air  were  utilized.  The  liquid  and  gas 
were  expanded  fro*  0  low  vsltx  .ty  and  a  prp-s^ure  of  500  peig  to  •ambient 
pt.-^wure.  Predicted  exit  velocities  ranged  between  ?0C  and  90°  ft/sec 
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for  the  air  and  between  600  and  800  ft/sac  for  the  water. 

Droplet  break-up  (or  critical  decomposition)  was  considered  in  the 
subject  analysis.  An  assumed  initial  droplet  diameter  of  0.C20  inches 
decreased  in  size  to  between  0,0008  and  0.0010  inches  at  the  nozzle  exit. 
The  Weber  number  governs  droplet  fcreak-up  and  a  critical  Weber  number 
of  6.3  was  employed  herein. 

Experimental  work  was  performed  in  order  to  check  the  theoretically 
predicted  pressure  profiles,  flow  rates  of  liquid  and  ga*»,  and  the 
fhrust.  The  theoretical  and  experimental  pressure  profiles  matched 
quite  closely  for  both  nozzles  investigated,  ffrisn  the  gas  flow  rate  was 
one-tenth  that  of  the  liquid  flow  rate  by  weight.,  the  pressure  profiles 
matched  more  closely  tlian  when  the  gas  flow  rate  was  0.15  that,  of  th< 
liquid  flow  rats.  The  predicted  total  flow  rates  were  within  four  to 
nine  percent  of  those  obtained  experimentally. 

The  thrust  was  predicted  to  an  accuracy  of  approximately  throe  per¬ 
cent.  for  the  short  nczxle,  but  to  an  accuracy  of  only  approximately 
twenty  percent,  for  the  long  noczle. 

From  the  investigation  it  was  concluded  that  tne  subject  analysis 
can  bo  utilised  to  predict  the  flow  characteristics  of  a  two-phase  flow 
of  liquid  drupe,  in  a  gas  stream  for  relatively  short,  males  with  the 
flow  r«lo  of  liquid  approximately  ten  times  the  flow  rate  of  gas  by 
weight. 


1  INTRODUCTION 


Die  analytical  method  developed  herein  can  be  employed  to  determine 
the  following  flow  chara-'terieties  of  a  two-phase  flow  of  liquid  drops 
in  a  gas  stream  as  the  mixture  expands  through  a  converging-diverging 
noizle-  gas  velocity,  droplet  velocity,  static  pressure,  and  droplet 
else.  Thir  --rk  was  initiated  in  conjunction  with  an  i  "vastigation  of 
the  operating  characteristics  of  a  gas-driven  jet  pump  (1)  (2)*.  The 
drive  notzle  for  such  a  device  ->erates  with  a  typically  low  mixture 
ratio.**  Previous  investigations  at  the  Jet  Propulsion  Center,  Purdue 
University  (1,2)  have  shown  that  an  operating  mixture  ratio  of  approxi¬ 
mately  0.10  should  be  utilised  for  optimum  jet  puap  operation.  At  lo*. 
mixture  ratios  a  larger  quantity  of  liquid  J.e  pumped  for  a  given  flow 
rate  of  gas  than  et  high  mixture  ratios. 

A  two-phase  nossle  is  a  device  through  which  two  fluids  are  accel¬ 
erated  froei  low  velocities  and  a  high  pressure  to  high  velocities  and  a 
low  pressure.  Ir.  the  subject  analysis  water  and  air  are  the  two  fluids. 
The  two  fluid*  are  injected  into  the  noatlc  from  an  injector  which  is 
designed  to  break  the  water  up  into  droplets  and  distribute  them  evenly 


*  Numbers  in  parenthesis  refer  to  references  appearing  in  the  rear  of 
the  report. 

**  Mixture  Ratio  -  defined  aethe  jaes  flow  rate  of  gas  (N^)  over  the 
meet  flow  rate  of  liquid  (W^)  -  or  v\ 
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in  tha  air.  When  the  two  fluids  enter  the  converging  portion  of  the 
nozzle  the  gas  expands,  and  therefore  its  velocity  increases  more  rap¬ 
idly  than  the  liquid.  The  liquid  droplets  accelerate  as  a  result  of 
(1)  the  decreasing  pressure  in  the  direction  of  flow  and  (2)  the  drag 
forces  exerted  on  them  by  the  faster  moving  gas. 

The  droplets  will  undergo  critical  decomposition* **  if  the  ir.  "tial 
forces  exerted  on  them  exceed  the  surface  tension  forces.  Thit, 
emon  is  governed  by  the  Weber  mmber.  The  Weber  number  is  'ined  as 
follows i 


N 


We 


(1) 


where 


»  gas  density 

Vjj®  -  square  of  the  relative  velocity  betwesn  gas 
and  liquid 
r  «  droplet  radius 
O'  *  surface  tension 

The  critical  Weber  number  .or  water  droplets  ir.  air  is  approxi¬ 
mately  6,3.  This  average  value  of  the  critical  Weber  lumber  wee 
determined  experimentally  by  Isshikl  for  water  droplets  with  diameter* 
tram  0.0181  in.  to  0.1924  in.^).  The  initial  drop  ditmetsr  assumed  in 
this  analysis  was  0.02  ih.  and  the  final  diameter  ranged  between  0.0008 
in.  and  0,0010  in. 


*  Critical  Dmcompoe ‘ 1  ion  -  Ekreaking  up  of  droplets  into  mealier  drop¬ 
lets. 

**  Weber  number  at  which  critical  decomposition  will  occur. 


lashiki's  work  indicated  that  the  critical  Weber  number  did  not  vary 
greatly  with  droplet  size.  No  other  data  could  be  found  for  smaller 
droplets  and  therefore,  the  average  value  of  6.3  determined  by  Isehiki 
was  utilized 

During  the  expansion  of  the  two-phase  mixture  through  the  nozzle, 
thermal  energy  is  transferred  between  the  liquid  and  the  gas.  In  the 
subject  investigation,  the  water  is  assumed  to  gi.c  up  heat  to  the  ex¬ 
panding  air.  The  amount  of  heat  transferred  between  the  phases  in  a 
given  period  of  time  depends  upon  (l)  the  difference  in  temperature  of 
the  two  phases,  (?)  the  amount  of  heat  lost  through  the  nozzle  boundaries, 
(3)  the  shape  and  size  of  the  liquio  droplets,  (U)  the  ’  locitiee  of  each 
of  the  fluids,  and  (5)  the  thermal  characteristics  of  the  individual 
phases  vfcich  includes  the  heat  transfer  coefficients  and  the  effects  of 
condensation  and  vaporization. 

There  ara  two  limiting  cases  of  the  hast  transfer  between  phases 
(l)  perfect  thermal  equilibrium  and  (2)  no  thermal  energy  transfer. 

Perfect  thermal  equilibrium  assumes  an  infinite  heat  transfer  rate  be¬ 
tween  phases  such  that  both  of  the  fluids  remain  at  the  same  temperature 
throughout  the  expansion.  No  thermal  aner.y  transfer  assumes  that  no 
heat  will  be  transferred  between  the  two  fluids.  In  this  case,  the 
temperature  of  the  gas  decreases  during  the  expansion  following  the 
lsentropic  r elationships  governing  perfect  gas  flow  through  a  nozzle. 

The  temperature  of  the  wtor  remains  constant  at  the  initial  inlet 
value. 

The  actual  case  of  heat  transfer  between  the  two  phases  lies  so»o- 
whsre  between  the  two  limiting  cases,  D.  G.  Elliott  and  D.  L.  Crabtree 


have  treated  theso  cases  in  detail  (i)  (2)r 

There  are  several  causes  Tor  the  loss  in  total  momentum  in  a  two- 
phase  expansion  Some  of  these  are  (1)  wall  friction,  (2)  dr*,  .lex  de¬ 
composition,  (3)  non-uniform  distribution  of  the  droplets  in  the  gas, 

(4)  non-ideal  heat,  transfer  between  the  two  phases  and  (3)  the  relative 
velocity  (or  slip)  between  phaseu . 

The  effxcts  of  the  above  mentioned  causes  for  a  i. ja  in  total  mo¬ 
mentum  of  the  mixture  have  been  investigated  in  previous  work  at  the  Jet 
Propulsion  Center,  Purdue  University  (l)  (2)  (4), 


| 
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2  TSSOflSnCAL  ARAU3ZS 

2-1  Purpose  of  Analysts 

The  subject  analysis  of  a  two-phase  expansion  of  liquid  drops  in  a  gas 
stream  flowing  through  a  converging-diverging  noss&Iva  ta  made  in  order  that 
the  droplet  velocity,  gas  velocity,  static  pressure  and  droplet  diameter 
may  be  predicted  a:  functions  of  axial  distance  along  the  nosrle.  More 
specif  lead,  ly,  the  subject  analysis  is  made  iu  order  to  determine  the  above 
mentioned  flov  parameters  in  a  mixture  containing  a  much  larger  amount  of 
liquid  than  gas  by  weight  (typically  *'r  »  10  Wr). 

2-2  Assumptlonc 

The  following  assumptions  were  employed  in  this  analysis : 

1.  The  tvo-pfc'vjc  flov  is  steady  end  one  dimensional . 

2.  Drag  forces  exist  between  the  liquid  droplets  and  the  g*fl. 

3.  There  arc  no  losses  das  to  wall  friction. 

4.  The  liquid  is  in  the  form  of  droplets  with  an  assumed  initial 
equivalent  spherical  diameter. 

5.  Critical  decomposition  of  the  droplets  (droplet  creak -up)  does 
occur.  The  droplets  are  uniformly  distributed  in  the  gke  uad 
U.-.S  of  equal  size  at  any  cross -section  ?  the  no ! tie . 

6.  Then*  is  no  interaction  between  droplets. 

7-  ’The  liquid  is  incompressible,  has  a  constant  specific  heat, 
and  has  no  vapo r  pressure. 
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6.  y*  gas  is  -srffect  and  hr*  cc ns bant  specific  h«at  at  constant 
pressure  and  constant  volume. 

,  So  external  work  i#  performed. 

10.  T'Oers  ia  no  beat  transfer  across  the  no* tic  boundaries. 

11.  Potential  energy  ia  a  constant . 

12.  She  flow  is  ieothcimi.. 

?|je  assumption  of  isotfearmal  flow  requires  further  explanation.  Tfcs 
equations  governing  tU®  tvo  limiting  cases  of  heat  transfer  between  the 
phases,  vhi.'-h  were  discursed  above,  have  bean  derived  by  ©.  0.  SXHott  Cl) 
Sb»  following  tv«>  equations  were  taken  fro©  tfcs*  uvek- 

Uniting  Case  I:  'fberaai  equilibrium 
between  phases  with  an  loa«  la  *»tn>Ry. 


r  «  x 
Tt  is 


Zn  i  t  ioll-.y .  *r 

1. 

Jtft«r  expansion,  T.  ■«  T„  *  t 


(w„/v„  kWjv) 

U  JjU 

c.  ♦  Cjw  /v  ) 

.7.  m  -’P  \  1  1 *' 

‘L  TG  /  \ 

XT'  *  JjT  "  i  l 
'S  ‘  S  ; 


wire 

f„  -  temperature  of  mixture  at  inlet  to  no tale, 

7,  »  liquid  temperature  at  any  cross  section  of  tbe 

i. 

T,  »  gas  twsperatuM  at  use?  «fo*s  section  of  the  nctsl*, 

G 

y  -  <*•»*  flow  rate  of  g»» . 

Q 

'«'  -  smlss  flow  rate  of  liquid, 

§  »  uBiversal  gas  core  it, at. 


(2) 
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W  *  molecular  weight  of  the  gas, 

J  »  mechanical  equivalent  of  heat,  778  ft  lb/3, 

CT  «  specific  heat  of  the  liquid, 

Cp  »  specific  heat  of  the  gas  et  constant  pressure, 

*  pressure  of  the  mixture  at  inlet  to  nozzle,  and 

P  >v  pressure  at  any  point  downstream  of  the  inlet. 

Limiting  Case  II •  No  heat  transfer  -between 

phase.?  and  no  I035  in  entropy. 

Initially,  T  =  I  -  T  . 

L  G  ;< 

Tr  *  c oust ant 

After  expansion,  T.  <  T. ,  11  *■  T„ 

\j  L  L  i » 


R/J  C>’ 
P  \ 


(3) 


The  inlot  pressure  and  temperature  of  t^s  two-phase  mixture  «tploy*d 

in  this  investigation  were  5  \  U .  7  psia  and  r',0°R  respectively.  The  drop  s 

in  temperature  which  would  occur  for  the  two  limiting  cases  are  found  by 

substituting  f.  >  51^.7  psia  and  T„  -  510°R  with  the  corresponding  specific 
n  n 

heats  and  gas  constant  into  equations  2  and  3  above. 


?  -  2 1* .  '  psia 

«  51L."1  pain 
. ' ' V  *  p3‘3  ft  lb/ lb  R 


for 
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CL  -  1.0  3/lb  R 
Cp  »  0.24  B/lb  R 

Ihe  rco^lts  are  shown  in  Table  1. 

Table  1 

Results  of  Calculations  for  Limiting 


Cases 

of  Heat 

Transfer 

Better;  Ifcas 

38 

Case 

“<A 

P 

T 

G 

T 

L 

*tg 

ATi, 

(psia) 

(R) 

I 

0  JO 

14.7 

496 

496 

12 

12 

J 

0.15 

14.7 

487 

487 

23 

23 

II 

any 

value 

14  7 

184 

510 

326 

Rxpeijjuental  vcrk  performed  at  the  Jet  PronulBion  Center,  Purdue 
University,  haa  indicated  that  the  real  case  of  heat  transfer  between 
phases  ie  closer  to  that  of  thermal  equllibriun  than  to  that  of  no  heat 
transfer  between  phases  (l)(2){4).  me  thermal  equilibrium  temperature 
drop  for  a  mixture  ratio  of  0.10  lr  only  12°R.  for  these  reaeons  the 
assumption  of  isothermal  flow  vas  made.  A  eo-ealled  "average  temperature" 
was  employed  throughout  this  analysis.  This  "average  temperature”  is  the 
average  between  the  inlet  temperature  chosen  and  the  exit  temperature  of 
the  mixture  calculated  with  the  thermal  equilibria*!  equation  2.  The 
thnrraal  equilibrium  temperature  drop  for  a  mixture  ratio  of  0.30  is  2^°R. 

It  should  therefore  be  expected  that  the  isothermal  model  will  hare  limited 


accuracy  us  the  aixt are  ratio  increase#  in  value  much  above  0.10. 


In  order  to  find  a  more  correct  drag  coefficient  for  the  droplet  than 
that  obtained  from  the  plot  of  Reynolds  Number  vs.  Drag  Coefficient  for  a 
solid  sixers  in  &  gas  stream  (5),  a  survey  was  made  of  existing  literature 
on  the  subject,  of  drag  on  droplets  in  an  air  stream.  It  was  decided  that 
the  best  data  to  utilize  would  be  that  experimentally  determined  by  E.  «»bin, 
R.  B.  Lawhead,  and  R.  Schallenmuller  at  Rocketdyne,  and  by  R.  D.  Ingebo 
at  the  La vis  Flight  Propulsion  Laboratory  (6) (7).  The  drag  coefficient  of 
liquid  droplets  having  diameters  from  0.004  in.  to  0.04  in.  was  determined 
to  be  approximately  one  and  appeared  to  be  independent  o~  Reynolds  number 
for  Reynolds  mssber  greater  than  fifty  (6).  The  Reynolds  number  referred 
to  throughout  this  work  is  the  droplet  Reynolds  number.  It  is  calculated 
by  employing  the  relative  velocity  between  phases  and  the  droplet  diameter 
in  the  defining  equation  for  the  Reynolds  number.*  For  very  small  droplets 
and  Reynolds  numbers  less  than  fifty,  Ingebo's  data  were  in  agreement  with 
that  obtained  at  Rocketdyne.  Ingebo's  equation  for  the  drag  coefficient 
of  a  small  droplet  in  an  air  stream  is  as  follows: 

£L_ 

0.84 


In  the  subject  analysis  the  drag  coefficient  was  assumed  to  be  one  fer 
Reynolds  nusbere  above  fifty.  For  Reynold"  Tiber  c  less  than  fifty, 

Ingebo’s  equation  for  the  drag  coefficient  .os  employed. 


*  See  Section  2-4.5 


Hie  following  equations  and  relationships  were  employed  in  the  subject 
analysis : 

1.  a  maaeutum  equation  for  the  system, 

2.  a  force  balence  on  the  droplet, 

3.  a  continuity  equation  for  the  mixture, 

4.  the  equation  of  3tste  for  the  gas 

5.  the  drag  equation, 

6.  the  defining  relationship  lor  Reynolds  number 

7.  the  defining  relationship  for  seber  mmber, 

8.  the  relationship  ol'  nozzle  flow  area  to  axial  distance  along 
the  nozzle,  i'nd 

9.  the  relationship  that  the  total  flow  area  is  equal  to  the  3um  of 
the  flow  areas  of  gas  and  liquid  at  any  cross -section, 

2-4  Derivation  of  the  Squat ion j  Smployed  in  the  Analysis 

2-4.1  Momentum  equaticJ|  for  the  system 

Consider  a  small  cross-section  of  the  nozzle  of  width  dx  shown  in 
fig.  1.  At  the  left  face  the  gas  velocity,  liquid  velocity,  pressure,  and 
area  are  denoted  by  V^,  V^,  p,  and  respectively.  At  the  right  face  all 
of  the  quantities  have  increased  by  differential  amounts  1  Indicated  in 
Fig.  1. 

Suasdng  the  forces  in  the  X-direction. 


x 


PAp  -  (p+^pMAp-HiO  *-  (?*dp/2)  sin a  (dAj/slna) 
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I  Multiplying  products  and  neglecting  differentials  of  higher  order  than  one, 

i 

i.  “  x  '"*T  dp  * 

The  sunmation  if  the  forces  in  the  X-direction  is  equal  to  the  change  in 
noosentun  of  the  g3S  and  liquid  in  thau  direction. 

2Fx  3-VP'¥VL+VVG 

Substituting  Hj.  =  and  K,  - 

I  rx  -yip  -  *  WA  w 

where 

*  mass  flow  rate  of  liquid, 

't,  *  mass  flow  rate  of  gas, 

0^  »  liquid  density, 
pQ  »  average  gas  dencity, 

Aq  -  average  flow  area  of  the  gas, 

A^  *  average  flow  area  of  the  liquid, 

»  average  liquid  velocity,  and 
«  average  g.  velocity. 

«  average  total  flcm  are* 

tv- arranging  Equation  4, 

.  Wl  *  Wowo 

dp  * - - -  (j) 

'  *T 

*  The  bar  (  )  over  any  quantity  indicates  the  average  value  of  that 

quantity  in  any  given  increment  4x. 
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2-4.2  Force  balance  on  the  droplet 

Consider  a  single  spherical  droplet  of  liquid  traveling  in  a  foster 
moving  expanding  gas.  Tvu  forces  tend  to  accelerate  the  droplet;  (l)  tK“ 
pressuro  force  resulting  from  the  pressure  gradient  in  the  direction  of 
flow,  and  (2/  the  drag  force  exerted  an  the  droplet  by  the  gas.  Writing 
a  force  balance  on  the  droplet, 

FP  ?D  "  “«°d  ’  VV"  =  Vl®!/"  (6) 

where 

Fp  «  force  due  to  preBsu.  j, 

?D  "  force  due  to  drag, 
a.  sui  of  the  droplet, 

U 

a.  *  acceleration  of  tb«  droplet,  and 
a 

VT  =•  droplet  velocity. 

The  farce  due  to  drag  is  determined  by  the  standard  drag  equation. 

h  m  (1/2)  W$i>  (7) 

where 

»  average  value  of  the  dr.vg  coefficient  for  the  it>crnmct  in 
question,  and 

A,  «  average  projected  area  of  the  drop.*t. 

U 

The  force  exerted  on  the  droplet  as  a  result,  of  the  preooure  variation 
around  it  can  be  determined  by  ccnaidering  a  small  portion  of  the  surface 
area  u“  the  droplet  4A  as  abovn  in  Fig.  2.  A  preaeure  force  P  act*  upon 
this  area  1A. 

The  droplet  hae  &  radius  r.  Tfcep,  dA  «  rdfi  (r  sin  8  d^)  and 
X  -  r  a  in  8  cov  0. 


u 


Tbs  pressure  variation  across  the  droplet  is  ass>^ed  to  be  linear  in  the 
direction  of  flow  (X-  action) .  Then, 


and 


P  «  (dp/dX)  X 


F  =  P  sin  9  cos  0 


The  pressure  fore*  acta  in  tae  negative  direction  in  the  quadrant  shown. 
Then, 


dJF  *  -  P  4A  *  •  •''P  sin  9  cos  0)  dA 

X  X 

*  -  (dp/dX) (X)  sic  9  cos  0  (rc  ain  0  d0  d0) 

»  -  (4p/dX)(r  sin  9  cos  0)  r2f'"‘‘0  cos  0  40  d0) 

=»  -  (dp/dX)  r3sin30  d0  cob20  d0  (8) 


Integrating  equation  8, 

Fx  -  -  (dp/dX)  r3^  ^  (sin30  40)  COB^  d0 


X  2« 


0  0 
x 


2% 


-  -  (dp/dX)  r3^  (sin30  d0)(0/2  sin 


2*/k) 


* 

(dp/dX)  rJ*  \  ain30  d0 
J0 

-  (dp/dX)  r  a  ^  cos  9  (sln20  +  c)~| 


*  -  (dp/dX)  r3*  (4/3) 

-  -  (voluee  of  droplet)  (dp/dX) 


therefore,  the  pressure  .  jrcs  in  the  X -direction  is, 


F  p  ■  (a p/ax) 


where 


dp/dX  is  negative,  and 

Vd  is  the  volume  of  a  droplet. 

Substituting  equations  7  and  9  into  equation  6  gives: 

Fp  +  Fd  =  -  (4/3)  *  r3  (dp/dX)  +  (l/2)  ^trr2^  =  md?L  (dVj/dX) 

Substituting  (h/3  n  r3^.  for  and  dividing  through  by  (4/3)  %  r3  yields, 
-  (dp/dx)  +  (3/8)  d^/F  «  (Jlvl  (dVL/dX)  (10) 


Rearranging  equation  10, 


dV.  *  zJlR  t  (3/8) 

*L\  W 


The  t  on  the  second  term  designates  whether  or  not  the  drag  on  the  droplet 
is  positive  or  negative. 


2-4. 3  Continuity  ec.ui tlon  for  toe  mixture 


WT  3  WL  +  WG  ”  *  lVl  4  S  gV7G 


W  *  ft  T  A.  V. 
L  T.  L  l 


UG  ■  Wo 


*  total  mass  rate  of  l'low, 

X-  =  specific  weight  of  tne  liquid,  end 

*  specific  weight  of  the  gas. 

2-4 . 4  Equation,  of  state  for  the  gas 

p  =  SG  RT  -  PqSRT 

where 

T  =  mixture  temperature  chosen,  and 
R  *  gas  constant  for  air. 


2-4.5  Defining  equation  for  Reynolds  number 


"a 


W* 


inertial  force 


viscous  force 


where 


d.,  -  average  droplet  diameter,  and 
u  =  viscosity  of  the  gas. 


2 


-4.b 


Defining  equation  for  Weber  cusaber 


N 


’We 


p..yDr  inertial  force 
u  n  . . 

<s'  surface  tens  ion 


Rearranging  equation  13, 

Ao  ' 


(15) 


(16) 


(1) 


*ad  differentiating  gives: 


"i  *  -  *  -  \  i*A> 
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Differentiating  equation  ifi  and  rearranging  terms, 


“o  -  “t  -  “l 


Differentiating  equation  15  with  T  held  constant  and  rearranging  terms, 


dpr  -  dp/ gRT 


Rearranging  equation  Ik,  differentiating,,  and  substituting 


yields 


jg  '  &cAA 


dVG  ^dp(/^G^  *  '~a 


Summarizing  the  above  equations: 


«A*t  4vt  *  Wo  "o 


4  A.VSCr.  ** 

dV  -  -  i  (  1/8)  - - 

*A  A? 


nh  -  WA 


-  \  (dvA} 


apjj  "  ap/RTg 


dVQ  -  Vq  (dpr/p0)  -  VG  (d^g/XJ 


(21) 


These  are  the  equations  in  the  foraa  utilises  to  determine  the  flow 
characteristics.  All  equations  ust  be  made  dimensionally  correct  oefcre 
calculations  are  performed.  These  equations  were  programmed  for  solution 
on  an  HPC  tOGO  digital  computer. 


2-5  Theoretical  Hozslea  Investigated 

The  gas  viscosity  and  liquid  surface  t*. us i or  were  considered  to  &a 
functions  of  tes^peratur*  only  and  ««  u4«.rsfor«  held  constant. 

Two  nozzle  ccafigaratioo*  v*»t  iovautigated  Sots  configure  5  ices  had 
the  seas  inlet,  throat  ,  u.nd  exit  areas,  throat  radius  zsA  cor.  'ergicj  argi*. 
I^ble  2  I'rts  tha  data  aegloped  ir,  tbr  analysis  and  the  physio  el  fcaensioos 
of  the  nozzle®. 

fig’.u'a  3  is  a  plot  of  total  aoale  flow  area  as  a  functicr.  of  a*i*l 
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AXIAL  VARIATION  OF 


2-6  Method  of  Solution 

Numerical  iteration  was  used  with  the  above  listed  equations  in  order 
to  determine  the  flow  parameters  as  functions  of  the  axial  length  of  the 
nozzle.  Small  increments  of  length  were  taken  along  the  nozzle  axis,  and 
through  each  increment  the  average  value  of  each  of  the  flow  parameters 
was  determined.  Hie  inlet  conditions  are  the  initial  conditions  for  the 
first  increment.  After  the  flow  parameters  have  been  determined  for  the 
first  increment,  the  values  of  the  parameters  at  the  exit  of  increment 
one  are  used  as  the  initial,  values  for  the  second  inert  'nt  and  so  on, 
until  the  exit  section  of  the  nozzle  is  reached. 

When  calculating  the  differentials  of  the  flow  parameters  with  the 
above  listed  equations,  overage  values  of  the  parameters  are  used  in  many 
of  the  expressions.  When  performing  the  ca.lcul.at ions  for  the  first  time 
for  a  given  increment,  the  average  value  of  each  of  the  parameters  is  not 
known  and  the r  fore  the  initial  v  ■ ues  of  the  parameters  ire  used  for  the 
average  values. 

After  the  differential  equations  have  been  solved  the  first  time  for 
a  given  increment,  an  average  value  of  each  parameter  can  be  calculated  an 
follows : 

Z  -  Z  *  dZ/2 

where 

Z  «  the  average  value  of  a  flow  parameter, 

*  initial  value  of  the  flow  parameter  at  the  beginning  of  any 
increment,  and 

dZ  «  dlf f-rential  value  cf  the  parameter  calculated  with  the  cor¬ 


responding  cq  .tion. 


The  equations  are  ag^.in  solved  for  the  same  increment,  hut  this  time 
using  the  average  values  vliit  have  been  calculated.  This  process  is  re¬ 
peated  aa  many  times  as  is  necessary  (for  each  increment)  to  ootaiu  the 
desired  accuracy  of  each  of  the  flow  parameters  at  the  exit  of  the  incre¬ 
ment.  The  value  of  a  parameter  at  the  exit  of  a  given  increment  is  found 
by  adding  the  differential  quantity  to  the  initial  value. 

Small  increments  were  taken  along  the  fixed  geometric  shapes  chosen 
and  therefore,  the  total  flow  area  for  each  increment  is  known. 

The  accuracy  to  which  the  parameters  are  calculated  effects  the  end 
results  since  ail  errors  are  multiplied  when  proceeding  from  increment  to 
increment.  All  calculations  in  this  analysis  were  made  with  an  accuracy 
of  0.’0£. 

The  size  of  the  increment  employed  greatly  effects  the  results  obtained. 
It  wa u  necessary  to  determine  the  smallest  increment  that  should  he  taken 
in  order  that  the  calculated  parameters  would  not  change  at  a  given  section 
of  "he  nozzle  If  a  still  smaller  rement  were  employed.  This  increment 
size  was  determined  to  be  approximately  twice  the  diameter  of  the  droplet 
in  length  at  any  position. 

Seven  significant  figures  of  accuracy  of  the  inlet  conditions  were 
not  sufficient  to  allow  the  expansion  to  proceed  tc  completion  at  14.7  psia. 
Tor  example,  let  the  correct  seventh  significant,  figure  of  one  of  the  para¬ 
meters  at  the  inlet  be  five.  If  a  trial  solution  is  attempted  using  four, 
the  solution  would  begin  to  divor^i  in  one  direction  (a  pressure  rise)  ct 
•cxse  point  along  the  nozzle  axis.  If  Instead  a  six  is  tried,  w.ve  solution 
would  d. verge  in  the  other  direction  (a  rapid  pressure  drop).  The  solu¬ 
tion  must  be  restarted  at  the  inlot  tc  the  nozzle,  with  the  new  chosen 


values  of  the  inlet  conditions  (VL,  VQ,  1^),  each  time  that  the  solu¬ 
tion  is  found  to  diverge  in  either  direction.  2he  RPC  4000  digital  com¬ 
puter,  when  employed  with  the  Purdue  Interpreter  routine,  could  only  he 
utilized  to  perform  accurate  calculations  to  seven  significant  figures. 

After  obtaining  the  correct  seventh  significant  figure  of  each  of  the 
inlet  conditions,  the  calculations  were  restarted  at  the  point  where  di¬ 
vergence  began  in  either  direction.  Integrated  forms  of  the  above  derived 
equations  v.'*e  then  utilized  in  an  iterative  form  as  discussed  fey  D.  L. 
Crabtree  (2). 

The  Weber  number  (equation  1)  was  utilized  to  determine  the  average 
droplet  size  in  any  particular  increment.  The  droplet  was  assumed  to  be 
of  constant  diameter  through  any  increment  at  the  average  value  determined. 
The  Weber  number  vas  calculated  at  the  nozzle  inlet.  If  this  calculated 
value  vus  less  than  the  critical  Weber  number  6.3,  the  droplet  remained 
conet ant  in  size  for  the  first  increment.  For  each  increment  the  Weber 
number  was  determined  and  if  it  remained  less  than  6.3  the  droplet  did  not 
undergo  critical  decomposition.  The  critical  value  is  the  maximum  value 
of  the  Weber  number  that  can  be  attained.  If  the  Weber  nuaber  becues 
equal  to  or  greater  than  6.3,  a  nev  droplet  radius  is  calculated  using  % 
value  of  6.3  for  the  Weber  number  in  equation  1. 

Figure  4  is  a  block  diagram  of  the  computer  pro,jrtw.  utilized  in  this 
analysis.  The  following  parameters  were  printed  out  for  each  increment; 

VL’  VQ>  r'  ^ye'  ?i«ure*  5,  6,  and  7  are  plots  of  the  predicated 

variations  in  the  flow  parameter*  with  axial  nozzle  length  Ibr  nozzles 
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3  EJGPIRriENTAT.  VERIFICATION 

3-1  Description  of  Experimental  Nozzles 

In  pnl«>r  to  verify  the  theoretical  predictions,  experimental  noz¬ 
zles  were  constructed.  The  experimental  nozzles  were  annular  in  shape 
and  mr«  constructed  with  the  same  area  vs,  axis  length  as  the  nosslus 
invast  igat  -ad  theoretic  ally. 

Figure  8  is  a  photograph  ox'  tha  component  parte  of  the  exiseriaental 
nozzle  ind  Figure  9  presents  a  cross-sectional  view  of  one  of  the  aa~ 
scnblatf  experimental  nozalas.  The  injector  was  the  same  type  as  uj„'  by 
1).  I,  Crabtree  \2) .  'lira  as.  ter  was  injected  through  small  tubes  with  so 
wiroler  flow  o-f  air  around  each  tube.  To  insure  that  the  strea,  ?  etait- 
tlng  froo  the  injector  ,«.  broker,  into  droplets,  a  stair, lass  steel  screen 
was  pis cad  in  i  he  path  of  flow  approximately  three  Inthas  from  the  in¬ 
jector  face.  Another  two  inches  were  allowed  before  the  converging 
portion  of  the  nozzle  began  in  order  to  stabilise  th<*  flow. 

>2  6s&d£U&jSLJl£  LiiSii^x 

Figure  10  1*  e  photograph  of  one  of  the  **oeri»*nUl  no tz lee  mounted 
on  the  thrust  stand,  The  1/S  in.  t..fc*s  which  a  -e  soldersd  to  the  outer 
nossie  casing  were  utilised  to  determine  the  pressure  profile  along  ttv- 
nossls  axii*.  Each  l/'<3  in,  tube  leads  to  a  0,060  in,  prsssodrs  tap  drilled 
through  the  outer  nossi*  casing,  Each  of  the  tubes  soldered  to  th<»  no- 
sis  casing  was  connected  to  a  common  manifold  through  resell*  valves. 
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VIEW  OF  THE  EXPERIMENTAL  NOZZLE 


s 


FIG  K)  EXPERIMENTAL  NOZZLE  MOULTED  ON  THRUST  STAND 
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The  pressure  at  a  particular  tap  was  measured  by  opening  the  correspond¬ 
ing  needle  valve  and  recording  the  manifold  pressure.  The  manifold 
pressure  was  measured  ty  a  Wlancko  pressure  transducer  and  was  recorded 
on  a  Brown  atrip-chart  recorder. 

Figure  11  ! s  a  schematic  diagram  of  the  facility  utilized  in  the 
experimental  investigation. 

The  nozzle  was  mounted  on  a  vertical  beam,  pivoted  at  the  upper  end 
as  shown  in  Figure  11.  The  thrust,  of  the  nozzle  was  measured  with  a 
Wlancko  force  transducer  and  recorded  on  a  frown  strip-chart  acorder* 

Water  was  supplied  to  the  nozzl6  at  pressures  between  5&3  and  1000 
psig.  from  four  tunka  which  have  a  5200  lb.  capacity.  The  water  flow 
rate  ws  measured  by  a  5/8  it  sharp-edgea  orifice,  The  pressure  drop 
across  the  orifice  was  met  *  red  with  a  Wiancao  pressure  transducer  and 
recorded  on  a  Brown  c'trip-chart  recorder.  A  milltameter  was  located  on 
the  control  panel  which  received  a  signal  from  the  Brown  recorder  to  In¬ 
dicate  the  liquid  i iow  rate.  The  liquid  flow  rate  was  controlled  with 
an  Annin  Domotor  valve  and  the  temperature  of  the  water  in  the  eupply 
Line  was  measured  with  an  iron-constantan  therooocouplo. 

Air  was  supplied  to  the  nozzle  from  several  banks  of  large  storage 
tanks  at  pressures  between  9^  and  1300  psig.  Grove  type  dose- loaded 
regulator  valves  w^-s  utilized  to  supply  the  nozzle  with  the  desired 
pressure.  The  ga*  flow  rate  was  measured  by  a  0.354  in.  sharp-edged 
orifice.  The  proesure  drop  across  ^he  orifice  and  the  pressure  upetraaa 
of  the  orifi.ee  were  measured  with  Wlancko  pressure  trnesducers  and  re¬ 
corded  on  Brown  strip-chart  recorders.  The  tessperature  of  the  air 
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upstream  cf  the  orifice  was  measured  with  an  iron-const-antan  thermocouple 
and  also  recorded  on  a  Brown  strip-chart  recorder.  Tha  air  flow  rate 
was  controlled  with  an  Annin  Domotor  valv-  and  was  indicated  on  the  flow 
panel  with  a  nilliameter  as  was  the  liquid  flew  rate. 

Both  liquid  and  gas  orifices  were  made  am'.  c«  .-.iterated  according  to 
A.S.M.E.  standards. 

The  pressure  at  the  entrance  of  the  converging  portion  of  t;.s  noz¬ 
zle  was  indicated  on  a  Bouraon-tube  pressure  gage  and  was  maintained  ct 
500  psig. 

The  pressure  and  thrust  transducers  were  calibrated  imedistely  be¬ 
fore.  and  after  each  experimental  run. 

3-3  Experimental  Data  Accumulation 

In  the  »xp ' rimental  invest igr •  on,  mixture  ratios  of  O.CS,  0.10, 
0.12,  and  0.1?  were  utilised  with  each  nozzle .  four  nozzles  were  uti¬ 
lized,  each  having  a  converging  angle  of  20°.  ”t verging  angles  of  7°, 
14°,  21°,  and  28°  were  employed.  Numerous  runs  were  made  in  order  to 
Insure*  that  the  results  of  the  particular  runs  desired  (those  identical 
to  the  theoretical  nozzle  calculations)  were  v; ) id ,  Data  taken  included 
the  foi.lowirg: 

1.  liquid  flow  rat*  , 

2.  ’»?  flow  rate, 

}.  thrust, 

i.  liquid  taripcrature , 

5.  gas  .empernture ,  and 

6.  pressure  r-rofil*,  along  nozzle  axis. 
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The  experimental  nozzles  which  were  identical  to  those  investigated 
j  theoretically  ar6  listed  in  Table  3* 


Table  } 

Nozzles  Employeu  in  Both 
Experimental  and  Theoretical  Analyses 


Nozzle 

Converging 

Angle 

Diverging 

/ngle 

Throat 

Arer. 

Inlet 

Area 

.Exit 

Area 

I 

20° 

7° 

0.267 

3-23/ 

2.450 

i: 

2C° 

21° 

0 .267 

3.237 

2.450 

i 


I 


i 


4  COMPARISON  OF  THEORETICAL 
AND  EXPERIMENTAL  RESULTS 


4-1  Pressure  Prefiles 

Figure  3  presents  the  flow  areas  of  the  two  nozzles  which  were  in¬ 
vestigated  both  experimentally  and  theoretically. 

Tn  Figures  12,  13,  end  14  the  predicted  pressure  profiles  are  com¬ 
pared  with  those  obtained  experimentally.  From  Figures  12,  13,  and  14 
it  can  be  seen  that  from  the  nozzle  inlet  into  the  diverging  section, 
the  measured  pressure  profiles  were  practically  identical  with  the  pre¬ 
dicted  profiles  in  all  three  cases  investigated.  In  each  of  the  three 
cases  the  experimental  pressure  profile  dropped  below  the  theoretical 
profile  in  the  diverging  portion  of  the  nozzle. 

In  the  diverging  portion  of  the  nozzle  both  phases  are  rapidly  ap¬ 
proaching  their  maximum  velocities.  In  this  region  of  high  liquid  and 
gas  velocities,  the  quantity  of  heat  transferred  in  a  small  finite 
period  of  time  is  minute  and  the  pressure  drops  faster  for  the  real 
case  than  for  the  isothermal  model.  As  the  pressure  drop'  per  unit 
length  becomes  small  in  the  latter  portion  of  the  diverging  section 
the  pressure  curves  again  match  quite  closely.  Ey  the  time  that  the 
mixture  has  traversed  the  entire  diverging  section,  the  liquid  has  had 
sufficient  time  to  transfer  some  heat  to  the  gas. 
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i-2  Flow  fiat a a 

A  comparison  of  the  predicted  and  measured  flow  rates  is  shown  in 
Table  4. 

Table  i 

Predicted  and  Measured  Flow  Rates 


Hozzle 

Mixture 

Ratio 

Predicted  To isl 
Flow  Rat.j  (lb/sec) 

Measured  Total 

Flow  Rate  (lb/sec) 

Percent 

Error 

y 

0.10 

10.04 

10.58 

5.4 

i 

0.15 

8.58 

9.35 

8.9 

ii 

0.10 

1C.  14 

10.58 

tO 

Tabi*  1  indicates  that,  the  isothermal  model  predicts  the  flow  rets 
mors  accurately  at  the  lower  mixture  ratios.  At  the  higher  mixture 
ratios  there  is  less  water  and  .more  air  snd  therefore,  less  heat  avail¬ 
able  for  a  larger  amount  ot  air  than  at  the  lower  mixture  ratios.  This 
increase  in  inaccuracy  of  the  isothermal  model  at  the  higher  mixture 
.-atioa  is  also  indicated  by  coaparir^  Figures  If!  and  13-  The  measured 
and  predicted  pressure  profiles  in  the-  diverging  portion  of  ftoisle  I 
match  more  closely  when  the  mixture  ratio  is  0.10  than  when  it  is  0.15. 
4-3  Thrust 

The  effective  n<mle  «xit  velocity  is  defined  as  follows: 

•  ♦ 

V„  *  Fg/Hr.  -  %  -  Pg/WT  •  Thrust /’Tot  vs;.  Flow  Rat- a  (22) 

and  the  a  viation  for  chruut  (F)  is  givon  by. 

F  -  1.1/,)  (V  •  Vo.’  *  <%  -  UU.  U3) 


.  K 

I- 
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idtere 

VGg  -  gad  exit  velocity 
“  liquid  exit  velocity 
p0  "  nozzle  exit  pressure 
Pa  -  ambient  pressure 
A„  ■  nozzle  exit  area 

O 

The  last  terra  of  Equation  23  can  be  dropped  since  p6  -  pa  in  both  ex¬ 
perimental  and  theoretical  cases.  The  predicted  and  measured  values 
of  the  effective  nozzle  exit  velocity  are  presented  in  Table  5. 

Table  5 

Predicted  and  !*asured  Values  of  the 
Effective  Nozzle  Qdt.  Velocity 


Nozzle 

Mixture 

Ratio 

Thrust 

(lb) 

Predicted 
Velocity (ft/ sec) 

Measured 

Velocity(ft/sec) 

Percent 

ESrror 

I 

0.10 

1 82.4 

699 

556 

25.7 

I 

0,15 

183.3 

778 

631 

23.3 

II 

0.10 

191*.  9 

620 

600 

3.3 

Froo  Table  5  it  is  seen  that  the  effective  nozzle  exit  velocity 
calculated  with  the  isothermal,  frictionless  model  was  far  in  error  for 
Houle  I.  The  calculation  for  Nozzle  II  however,  was  quite  accurate. 
Nozzle  I  had  a  diverging  angle  of  7°  and  an  overall  length  of  8.85  in. 
and  Nozzle  II  had  a  diverging  angle  of  21°  and  an  overall  length  of 
5.1  in.  The  results  presented  In  Table  5  indicate  that  friction  clays 
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a  dominate  roll  in  the  determination  of  the  gas  and  liquid  velocities 
when  long  nozzles  are  being  investigated.  Friction  between  the  phases 
and  between  the  miy*  .  e  and  th6  nozzl6  boundaries  was  neglected  in  this 
analysis.  Non-ideal  thermal  equilibrium  also  causes  a  reduction  in 
thrust.  (2) 


I 
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COLLUSIONS 


The  use  of  th®  frlcUoniuss,  isothermal  model  appear*  to  to  limited 
tc  those  type  nozzles  suited  for  ,Ut  pump  operation#  Those  non*!**  are 
short  to  conserve  weight,  and  operate  witn  '?#  mixture  ratios  in  order 
that  ••  much  liquid  be  pumped  as  posaiois  with  a  corresponding  low  ?tes 
consumption. 

Two  modifications  should  be  .made  to  the  isothermal  modal  in  order 
to  widen  the  -eng*  of  its  usefulness;  (1)  consideration  should  be  given 
to  the  heat  transfer  that  occurs  between  the  liquid  snd  the  gas.  and 
(2)  the  effects  of  wall  friction  should  be  included, 

An  analysis  which  takes  into  account  the  heat  transfer  between  the 
liquid  and  the  g» a  would  employ  the  following  equetions  and  relation¬ 
ships  in  addition  to  thoss  employed  in  the  ieotheraal  model. 

1.  Energy  equation  for  mixture 


fc^d?  ♦  V  (dV  / 

Jg>|  *  *  Id h  «  V  { 

dV  /Jg) 

IP  G  G  0 

J  L  L L  t 

i. 

2.  Tl  -  tihj 

3.  Heat  Balance  on  the  Droplat 

„i/dt  -  -\tJC  •/‘aYt  diVdt 

’  /W,  ih^a 

U.  &p*rla>anta  1  plot  oi  Nussalt  number  vs.  Reynolds  number 
for  spherical  droplets  (8). 

I 


where 


“  droplet  volu&a 

Ad  ■  droplet  surface  area 

h„  *  convectiva  he  transfer  coefficient 
c 

-  enthalpy  of  the  droplet 

q^j  »  convective  heat  transfer  rote  "rorr,  the  droplet 
AT  *  temperature  difference  between  the  gas  and  the  liquid 
•  droplet  density 

Equation  B-5,  which  is  derived  in  Appendix  B,  shows  the  effect  of  a 
temperature  drop  of  the  gas  on  the  gas  velocity  for  t'»a  adiabatio  model. 


dV  -  dVr 

u  adiabatic  u  isothermal 


model  model 

dig  is  negative  and  therefore  dV^  is  'ess  for  a  given  increment  when 
employing  the  adiabatic  model  than  wh,,i  employing  the  isothermal  model. 
A  lower  exit  gas  veloclt;  results  when  considering  the  flow  to  be  adi¬ 
abatic  which  in  turn  reduces  tho  thrust  that  can  be  obtained  from  the 
notice. 

5-2  m£U2n^j?xi;^«ai 

It  appears  that  o  friction,  coefficient  must  be  experimen*' ally  de¬ 
termined  for  the  type  of  two- phase  flow  discussed  heroin.  There  is 
friction  between  the  gas  and  liquid  and  between  the  mixture  and  the 
not*, l»  boundaries.  For  each  different  git.  und/or  liquid  used,  the  in¬ 
fluent.  of  friction  would  bo  different.  <*  theoretical  detamination  of 
the  friction  coefficient  appears  to  be  Improbable . 


5-3  fagwttfrti^gat 

The  analysis  presented  herein  should  be  programmed  for  solution  on 
a  fsater  digits!  computer.  If  the  hast  transfer  betw»«n  the  liquid  and 


ib 


t.he  gas  is  to  be  considered  it  would  be  imperative  that  „  faster  computer 
be  utilized.  To  determine  the  inlet  conditions  to  seven  signific  t  fig¬ 
ures  for  one  nozzle  and  mixture  ratio,  approximately  sixty  hours  of  com¬ 
puter  time  were  necessary  when  employing  the  RPC  4000,  The  solution 
would  take  a  considerable  amount  of  time  even  if  a  faster  computer  were 
utilized.  Another  method  of  numerical  solution  of  the  oquations  should 
therefore  be  investigated. 
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APPENDICES 


AK'PEMDIi  A 


NOJOCLATilW 

4d  “  acceleration  of  droplet,  ft/sec® 

Ajj  *  projected  area  of  droplet,  in.2 

Afl  »  noizle  exit  area,  in.2 

Aq  •  flow  area  of  gas,  in,® 

;W  -  flow  area  of  liquid,  in.~ 

*>  ono« dimensional  surface  area  of  diverging  nozA.lt,  in. 

Aq.  «  total  flow  area,  In.8 

Cp  “  dreg  coefficient 

Cp  *  fic  heat  of  water,  BAo-°il 

Cp  ■  specific  heat  of  air  at  constant  pressure,  B/lb~°R 

T  •  thrust t  lb 

Fp  *  drag  force,  lb 

F p  *  pressure  force,  lb 

¥  *  force  in  Indirection,  lb 

x 

g  m  gravity  constant,  32.174  ft/ sec* 

•J  -  technical  equivalent  of  heat  ,  778  ft-lb/B 
of  droplet,  slugs 

• 

4^  ”  ffiflss  flo.»  rate  el  gas,  aluga./sac 

a 

**  «a#s  flow  -at*  of  liquid,  siugs/aec 
Mjj  •»  ^aynold*  ruaber 
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%e  a  numoer 

p  “  pressure,  paia 

P  »  pressure,  psia 

p_  *  ambient:  pressure,  psia 
p  m  notzi  exit  pressure,  paia 
Fjj  *  noazle  inlet  pressure,  psia 

P  *  X-ccnponent  of  pressure  acting  on  droplet,  psia 

X 

ft  “  Universal  gas  const*!-.* ,  ft~lb/  stole -R 

R  »  gas  constant  for  air,  51*3  ft-lb/lb-R 

r  *  droplet  radius,  in. 

Tq  «  gas  temperature,  R 

Tj  23  liquid  temperature,  R 

Tjj  ■  inlet  mixture  temperature ,  ft 

Vq  "  gas  velocity,  ft/sec 

VL  »  liquid  velocity,  ft/sec 

V„  ■  effective  noxslo  exit  velocity,  ft/cec 

VR  -  relative  velocity  between  gas  and  liquid,  ft/sec 

W  •  molecular  weight  of  gas 
• 

Wg  -  mss  flow  rate  cf  gas,  lb/aoc 
■  i&ass  flow  rate  uf  liquid,  lb/ sec 

• 

“  total  flow  rate,  lb/sec 

1  ■  diatance  along  nozxle  axis,  in. 

2  “  any  flow  paraaetsr 

*  specific  weight  of  liquid,  lb/ ft* 

-  specific  weight  of  ^s,  lb/ft® 
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■  T 


M  *  gas  viscosity,  slugs/ft-otc 
»  density  of  liquid,  slug? 

^<0  "  density  of  gas,  slugs/ft8 
*y-H  «  droplet  volu&e,  in*3 
CT  **  surface  tension  of  liquid,  lb/ ft 
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AFPSfiDIS  B 


CO.' tF  .ARISON  OF  .4Di;,6,vTIC 


and  isotherual  equations 


In  order  to  compare  the  equations  for  the  adiabatic  and  the  iso¬ 


thermal  mi.de la ,  ail  of  the  equations  for  both  cases  are  left  alike 


except  for  the  equations  for  riV^. 


Rewriting  aquation  21 


VG  isotoeraal  '  (dV'V  “  %  G) 


Writing  the  "...  'vet  gss  law  for  the  adiabatic  model  (T  /  constant)  and 

G 


ra*rrar4i.ng  terms 


/Og  -  p/gRT 


Performing  1  of arithmic  difforention, 


W/>G  r  d-r/P  -  dVTG 


Substituting  equation  {&-•!)  into  equation  (21) 


udiabatic 

model 


'W3^  -  ^(dp/^dVkJ-' 


For  isothermal  .flow  cST„/T..  •  0  and  freo  equation  (B-l) 


d,%>G  '  W* 


Substituting  equation  5-3  into  equation  21 


dV,  -  -?«(dWX*)  "%  t<Wp) 

o  ieothansol  -  u 


V’  vy 


1  *  V  *V 


*  ft 


-if 


I. 

A. 

hs  • 


' 


r  ■■■«' 
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Substituting  equation  B-4  into  equation  B-2 


G  adiabatic 
model 


dV  *  V  (dTr/Tj 

G  isothermal  G  G  a 


4  f| 


-  *■  't  i- 


M. 
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"Diatribe*"1  of  this  report  iu»»  bean  aade  in  accordance  *ith  tht  Joint 
Anay-Ravy-Air  Korea  liquid  Propellant  Mailir^  Liat  of  March  1962.” 
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